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Abstract 

The scaling properties of correlation functions of non-scalar fields (con- 
structed from velocity derivatives) in isotropic hydrodynamic turbulence 
are characterized by a set of universal exponents. It is explained that these 
exponents also characterize the rate of decay of the effects of anisotropic 
forcing in developed turbulence. This set has never been measured in 
either numerical or laboratory experiments. These exponents are impor- 
tant for the general theory of turbulence, but also for modeling anisotropic 
flows. We propose in this letter how to measure these exponents using 
existing data bases of direct numerical simulations and by designing new 
laboratory experiments. 

Fundamental studies of turbulence tend to stress the model of isotropic, ho- 
mogeneous turbulence, and most theories and experiments since Kolmogorov's 
seminal work of 1941 |lj considered the universal (anomalous) exponents that 
characterize the isotropic characteristics of turbulent flows (see for example [El [| 
for recent reviews). In fact, most turbulent flows are not forced isotropically, 
and moreover even in isotropic flows there are important fields that are con- 
structed from velocity derivatives that transform under rotation as vectors or 
tensors rather than scalars. It has been known for quite a while that the second- 
order structure function (that depends on one separation vector) becomes more 
and more istotropic as the the separation scale goes down (see below). More- 
over, the rate of this isotropization process is governed by a universal exponent 
[H II [| • I n recen t papers [|[ ^) it was pointed out that this exponent is one 
of an infinite family of universal anomalous scaling exponents that were never 
considered in experiments and numerical simulations. Moreover, it was shown 
that in the context of passive scalar convection the anomalous scaling expo- 
nents that characterize the scaling properties of anisotropic fields also govern 
the rate of isotropization of the properties of the flow in the cascade process 
down to smaller and smaller scales [pi. In this Letter we show that the same 
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connection between the exponents of tensor fields and the exponents governing 
the rate of isotropization exists also in Navier-Stokes turbulence. This demon- 
stration follows from the assumption of "weak universality" of hydrodynamic 
turbulence that is discussed in jlO| . Theoretically the exponents discussed here 
are related to the appearance of the inner, viscous scale rj as a renormaliza- 
tion scale in addition to the more commonly known appearance of the outer, 
integral scale of turbulence L [jll| [l^, |J. In addition to being of fundamental 
interest these universal properties are also of importance in modeling realistic 
flows which are not isotropic. In this Letter we also propose how to measure 
these universal exponents in numerical and in physical experiments. We will 
focus here on stationary, space homogeneous turbulence which is however not 
necessarily isotropic. 

The simplest statistical quantity that is built from the fundamental velocity 
field u(r, t) that displays important contributions from anisotropy is the second- 
order structure functions of velocity differences w(ro|r,i) = u(r,£) — u(ro,t): 

S 2 (R) = (|w(r |r,i)| 2 ) , R = r — r . (1) 

where (...) stands for a suitably defined ensemble average. Due to space and 
time homogeneity 5*2 is a time independent function of the vector R. In isotropic 
turbulence the scaling properties of S2 (R) were widely discussed ^ ||, |l3| 

£ 2 (i?)-( £ -i?) 2/3 (fy , 

where e is the mean energy flux per unit time per unit mass, and 5 is the devi- 
ation of the scaling exponent C2 of the structure function from the Kolmogorov 
1941 (K41) prediction £2 = 2/3 — 6. In anisotropic turbulence 5 2 (R) depends on 
the direction of R, and we can decompose it into spherical harmonics according 
to the "multipole expansion" 

00 

5 2 (R) - ]TS 2 /(R), (3) 

;=o 

%(r) = J2 / s 2 (Ri)Y lm (i)di , (4) 
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where £ is a unit vector. In a scale invariant situation every component S2J 
scales like 

2 / 3 fR\ e 



S V (R) ~ (my lf\ cx R* , (5) 

where Si = f3i — 2/3. Comparing with Eq.(|^) we recognize that in this notation 
8 = 8q. The full spectrum of exponents Pi was found analytically || Q in the 



context of Kraichnan's model of passive scalar convection) 15 . For Navier-Stokes 
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turbulence (3 2 can be computed using perturbation theory ||, |?J (which disre- 
gards the nonperturbative effects leading to anomalous scaling |T^]) with the 
result f3 2 = 4/3. The corresponding result for /?o = C2 is 2/3, which experimen- 
tally is known to be of the order of 0.7 |l6|, The theory indicates that such 
deviations from the naive predictions stem from non-perturbative effects. It is 
likely that the perturbative result for the exponent [3 2 holds to a similar accu- 
racy. The large difference between (3 2 and (3q explains why isotropic scaling may 
be observed in anisotropic experiments; the contribution of S 2 , 2 to S 2 peels off 
rather quickly when R <C L. We do not possess any numerical estimates for the 
higher order values of [3i , but we expect them to be all positive and increasing 
with I. 

In the context of passive scalar convection we demonstated that the very 
same exponents f3e have an important role in the context of isotropic turbu- 
lence when we considered statistical quantities that depend on more than two 
coordinates]^. We point out here that the same is true for Navier-Stokes turbu- 
lence. Consider for example the correlation function of four velocity differences 

5 4 (Ri,R 2 ) ^(IwtrolrOlVMr^l 2 ) , (6) 

where Ri = ri — ro and R2 — r 2 — tq. As usual we assume that this, and 
all other correlators, are scale invariant. Mathematically this means that they 
are all homogeneous functions of their arguments as long as these are in the 
"inertial range". In other words S^ARi, AR2) = A^S^Ri, R2) where £4 is the 
scaling exponents of the 4'th order structure function: ^ |w(r |r!) | 4 ) oc iq 4 . In 
isotropic turbulence S^Ri, R 2 ) depends on the separations Rx, R 2 and on the 
angle 8 12 between these two vectors. We are interested in the limit R\ <C R 2 , 
but i?i and R 2 are both in the inertial interval. It was shown in that in 

this limit the leading dependence of S*4(Ri,R2) on R\ and R 2 is independent 
of 9i 2 and that it scales like 



Ri 



C2 



S 4 (Ri,R 2 ) oc l^j i# , Ri « R 2 . (7) 

In order the extract the sub-leading dependence on Ri we use a mulitpole 
decomposition of S4 in a way similar to (Q): 

00 

54(Ri,R2) = ^54^(Ri,R2), (8) 

1=0 

S 4 ,<(Ri,R 2 ) = J2 YURt) J S 2 (R 1 ili 2 )Y lm {i)di . 
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The generalization of (Q) to the anisotropic contributions is 

RxV* 



S v (Rl,R 2 ) oc ( j^-j Rf, R\ « R 2 . (9) 
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Note that there are two statements made here. The first is that the over-all 
exponent for this quantity is £4- This directly follows from the property of scale 
invariance. The second statement is that the scaling exponents characterizing 
the i?i-dependence of are the scaling exponents of S2/. This result follows 
from the assumption of "weak universality" used in the derivation of the fusion 
rules in |To| . Physically it is equivalent to the statement that the measured 
scaling exponents of the structure functions in turbulence are independent of 
the precise driving mecahnism at the outer scale of turbulence L. The existence 
of i?2-scale eddies and their effect on the statistics of the much smaller Ri 
eddies is similar to the existence of L-scale eddies and their effect on the scaling 
exponents in the inertial range. The velocity difference measured at points 
ri and ro is effected by the large eddies characterizing the velocity difference 
across R2 in a way that is similar to the effect of the boundary conditions 
at L on the velocity difference measured in 52- Similarly the exponents (3e 
characterize the rate of isotropization of S4 as a function of Ri. We stress 
that this isotropization is relevant in isotropic turbulence since we have in this 
function a built-in direction R2. When i?i is of the order of R2 the dependence 
on the angle 6*1,2 is all important. When Ri decreases this dependence weakens 
at a rated determined by the exponents Pi — /?o- 

Next in order of complication we consider S^Ri, R2, Ro) defined as 

5 4 (Ri,R 2 ,Ro) = (|w(r |n)| 2 |w(r |r2)r) , (10) 

where Ro ~ r' Q — Tq. This is a function of three separation and the three angles 
#i,0> #2.0 an d 0i,2- As before represent this function as a double multipolc- 
expansion with respect to the directions of Ri and R 2 : 

54(Rl, R2, Ro) = ^ S , 4,£i,£ 2 (RljR2 7Ro) • (11) 

In the limit R\ , R2 *C Ro these functions exhibit a universal scaling form similar 
to® 



S4A/ 2 (Ri,R2,Ro) oc ( J h£) ^ ■ (12) 

Finally we discuss correlations of anisotropic local fields constructed from 
velocity derivatives. These can be obtained by a limiting procedure starting 
from the fusion of two points in the vicinity of r as in (^) or in the vicinity 
of ro and r' Q as in (p^). As shown in M the simplest representatives of such 
fields contain a product of two derivatives daUpd^ug. Higher tensorial fields are 
obtained by taking additional derivatives from this field. To get clean scaling be- 
haviour we need to decompose these fields to combinations that give irreducible 
representations of the rotations and inversion group 0(3). Every irreducible 
representation is characterized by an index £, has a dimension (2£+ 1), and the 
(21 + 1) fields that form its basis transform like the spherical harmonic Yi :in . 
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The low orders representations are constructed with the help of the strain ten- 
sor Sq/3 = [dua/drp + duf3/dr a \/2 and the vorticity field u> a — ea^dup/dr^ 
(where e Q /3 7 is the fully antisymmetric pseudo-tensor). There are two scalar 
fields, Og 1 ^ = uj a Lu a and = s 2 — s a psp a each of which is a basis for 

one-dimensional irreducible representation with £ = 0. The pseudo-vector 
Of = s a pujf3 is a three-dimensional basis for an irreducible representation with 
1=1. There exist three traceless tensor fields each of which is a five-dimensional 
basis belonging to £ — 2 and taking care of 3x5 = 15 components. An example 
is 

0^ Q/3 (r) = SQ7 (r)s 7/3 (r) - S a0S 2 (r)/3 . (13) 

In addition we have one 3-rank pseudo tensor corresponding to £ = 3 and one 
4-rank tensor corresponding to £ = 4. The last two fields exhaust the remaining 
7 + 9 components We note that all the tensor fields Og are dimension- 
ally identical. However, dimensional analysis misses the point, and fields that 
transform differently under the symmetry group have different scaling expo- 
nents. The correlation functions (0^(r + R)O^(r)) all have different scaling 
exponents depending on £ and £': 

K U ,(R) = ((^(r + RjO^Cr)} oc iT"«' , (14) 

independent of n and n' . In particular the prediction is that the three cor- 
relation functions involving the scalar fields s 2 and to 1 have the same scaling 
exponents known as the "intermittency exponent" /i. Similarly the six correla- 
tions involving O2 have the same exponent (different from Note that the 
rank of K«/ is the sum of the ranks of the tensors in the correlation. 

At this point we want to discuss how to set up possible experiments to 
measure the new universal exponents f3g. Given a direct numerical simulation 
with anisotropic forcing, the most straightforward way it to simply compute 
S2,e(R) from the definitions ([!]) and ([|) and then to plot log-log plots of S-z,e 
vs. R, or even better, following the ideas JIt]] of "extended self-similarity" , 
of S2/ vs. /f>2,£=o- It is impossible to follow this route in standard laboratory 
experiments since the detailed angular information is not usually available. One 
can estimate the exponent P2 in anisotropic flows by measuring for example the 
longitudinal and transverse components of the second order structure function, 
and form a combination that vanishes in isotropic flows. Such a combination 
scales with R and the leading contribution is -fir 2 . This type of measurement 
was performed, see for example jll| and discussed in detail by Nelkin ||. The 
experimental evidence is that the numerical value of fa is indeed rather close 
to 4/3. Our point in this Letter is that the very same exponents @i play an 
important role also in isotropic flows. 

In laboratory (and also in atmospheric) experiments it is difficult to resolve 
the dissipative scales, and the direct measurements of the local fields Oi(r) is 
quite hard. Accordingly we will suggest experiments that are based on finite 
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differences in the incrtial range instead of gradients. Consider an experiment 
with a mean flow (like a wind tunnel or an atmospheric boundary layer) . Assign 
the direction of the mean flow to the rr-coordinatc. The minimal experimental 
set up calls for two local probes (like hot wires) positioned at r = (0, 0, 0) and 
i"i = (0, A, 0), separated by a distance A in the y-direction which is orthogonal 
to the mean flow. Under the standard Taylor hypothesis differences in time are 
interpreted as differences along the longitudinal x-direction. This means that 
one can measure the longitudinal projections au x (x, 0, 0) and bu x (x, A, 0). The 
coefficients a and b were introduced in recognition of the fact that in realistic 
experiments the two probes cannot be perfectly calibrated. Define now the 
longitudinal and transverse velocity differences 

w\\{x,A) = u x (x + A,0,0)~u x (x,0,0), (15) 
w ± (x,A) = u x {x,A,0) -u x (x,0,0) . (16) 

Next one can measure the longitudinal and transverse structure functions for 
A-separations 

S 2 \\(A) = (tf\(x,A)), S 2± (A) = ( W 2 ± (x,A)) . (17) 

In isotropic conditions these two quantities are related |D| by 

S 2 ±(A) = S 2 ||(A) + AdS 2 ||(A)/2dA , (18) 

and one can use this relation to assess the degree of isotropy on the scale A. 
Next we introduce the normalized squared of velocity differences in which the 
calibration constants are eliminated: 

w?,(x,A) wf,(x,A) 

"fi^-fesp w ^ s iktw (19) 

Finally we define two fields ^± according to 



¥+(x,A) = W]_(x, A) + Wn(x, A) 



*_(x,A) = Wl(x,A) - W { 2 (x,A) . (20) 

These fields have zero mean by construction, and we propose to measure their 

correlations K ++ (R), K--(R) and (R) across a scale R such that rj < A < 

R < L. The theoretical prediction is that the leading scaling form is 

K ++ {R) = (y + {x + R,A)t>+(x,A)) oc R->*++ , (21) 
K + _{R) = (*+(a; + i?,A)*_(a;,A)> oc , (22) 

K—(R) = (^^(x + R,A)^-(x,A)) ocR-v— , (23) 

where the exponents are 

^++ = 2/3 - C4 i Po = C2 , 

= Po + (h - (i , H— = 2/3 2 - C 4 • (24) 
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The reason for this prediction is that the field \&+ has a large projection on the 
zero'th spherical harmonic !b,m, whereas the field , 5- has no such projection 
by construction. Moreover the field also has no projection on li, m . It 
does have projections on 52, m and higher order spherical harmonics. Thus the 
leading scaling exponents appearing in correlations of is 02, whereas the 
leading scaling exponent appearing in correlations of ^ + is /3q, and Eqs.fcl 



23J) follow directly from (pL2|) . Note that the correlations in (21)-fl23|) are all 
dimensionally identical; yet we predict very different scaling exponents. This 
is just another way to explore the breakdown of dimensional analysis in fully 
developed turbulence. Using known experimental data (l(],[l7) (2 — 0.7, £4 — 1-2 
and our guess that (3% ~ 4/3 we expect: 

,i++~0.2, M __~1.4-1.5, H+- = /X+ + + /X — , (25) 

The last relation is asymptotically (in Re) exact. Note however that for a 
finite extent of the inertial interval sub-leading scaling contributions may be 
important and have to be carefully assessed. Nevertheless, the wide disparity 
between these scaling exponents promises a worthwhile experiment even if the 
inertial range is of the order of one decade. 

Direct numerical simulations offer additional ways to examine the correlation 
function K«<(R) ( |l4| ) of the anisotropic local fields. It was pointed to us by 



R. Benzi [19 that in numerical simulations the correlation function of gradient 
fields have a scaling range that is too short. Indeed, even the correlation function 
of the scalar dissipation e(r) = 2vs 2 {r) is not readily available. It was suggested 
that a better scaling behaviour is exhibited by integrals of the dissipation field 
over balls of radius R. The calculation is achieved by first considering 



£ * (ro) = 4^ / dr £(r) ' (26) 

from which one computes the time and space average 

(4) = ([ efl ( ro )-ef) • (27) 

For R in the inertial interval 

(e\) = C„K ee (R) oc iT" , (28) 

where C M ~ 1 and K ee (R) = ([e(r + R) — e][e(r ) — e]). Following the same 
idea |nj instead of computing in simulations the correlations K«/(R) we can 
measure the mean value of the product of two i?-ball integrated fields 

O^(r ) = / dr O t (r) . (29) 

AttK* J\r-r \<R 
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Accordingly 

<O^(r )O^(r )) = C«>K«,(R) cx iT"«' . (30) 

It should be stressed here that we did not offer yet any numerical estimates 
for /j,££i . The reason is deep: the local fields pick up a dissipative scale when 
the gradient is computed. In Navier-Stokes turbulence there is a multiplicity 
of dissipative scales that have non-trivial dependencies on the inertial-range 
separation distances which the relevant correlation function depends on. The 
full discussion of this issue will be available elsewhere jH^. Here we will just 
present theoretical prediction without derivation. They are: 

m . =2-(e + f3 i + p i , -2p . (31) 

For £ = £' — we recover the well known "bridge relation" which followed from 
the Kolmogorov refined similarity hypothesis, /j, = 2 — £g. The predictions for 
non-zero £ are novel and await confirmation. We stress that such a measurement 
can give information about /3i which is associated with the pseudo-vector field 
SafjWp. This exponent is not available from the rate of isotropization of S^R). 
It can be seen if the flow field does not have inversion symmetry but one needs 
to form a nonsymmetric second-order correlation function like 

(u a (r + R)u p (r)) - (u a (r - R)u/j(r)> oc R h . (32) 

Since this object is manifestly odd in R it vanishes when there exists inversion 
symmetry. Otherwise its leading scaling exponent is j3\ . This exponent is related 
to the existence of the flux of helicity and standard K41 arguments lead to the 
value 0i = 1, see for example This holds probably to the same accuracy 
as other K41 arguments. The correlation functions discussed above allow the 
measurement of this exponent in the presence of inversion symmetry. Even 
thought K«/ will vanish for £ even and £' odd, the correlation Kn is non-zero 
in any case. 

In summary we presented briefly the ideas that related the infinite set of uni- 
versal exponents characterizing the rate of isotropization of the simple second- 
order structure function under non-isotropic forcing with the same set of expo- 
nents that determines the scaling behaviour of correlation functions of tensorial 
anisotropic local fields in isotropic turbulence. The central role that these ex- 
ponents play warrants their measurement in laboratory and numerical experi- 
ments. We thus offered some simple ways to measure the low order exponents 
Pi and @2 m realistic experiments. We presented an estimate of the numerical 
values of these two exponents. The calculation of these exponents from first 
principles is a different task that is outside the scope of this Letter. 
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